1 

Correlation between Cirrus Particle Optical Properties: 
Microphysics and Implications for Spaceborne Remote Sensing 

J. Reichardt 1,2 , M. Hess 3 , S. Reichardt 1 ’ 2 , A. Behrendt 4 , T. J. McGee 2 

1. Joint Center for Earth Systems Technology, University of Maryland 
Baltimore County, Baltimore, MD 

2. Atmospheric Chemistry and Dynamics Branch, Laboratory for 
Atmospheres, NASA Goddard Space Flight Center, Greenbelt, MD 

3. Deutsches Fernerkundungsdatenzentrum, Deutsches Zentrum fur Luft- 
und Raumfahrt Oberpfaffenhofen, WeBling, Germany 

4. Radio Science Center for Space and Atmosphere, Kyoto University, Uji, 
Kyoto, Japan 


Submitted to 
Science 


Corresponding atithor address: Dr. Jens Reichardt, Code 916, NASA 
Goddard Space Flight Center, Greenbelt, MD 20771; Phone: 301-614-5989; 
Fax: 301-614-5903; Email: reichardt@code916.gsfc.nasa.gov. 


3 


■r 


Cirrus clouds are widely recognized as a major component of the energy 
budget of the Earth-atmosphere system ( 1 ), yet the magnitude of cirrus 
climate forcing remains highly uncertain, despite considerable efforts to 
quantify it ( 2 , 3 ). The difficulties stem primarily from the fact that the 
radiative forcing of cirrus clouds, which are predominantly comprised of ice 
particles, is crucially dependent on the particle shapes, sizes and number 
concentration. An accurate assessment of cirrus cloud forcing poses a m.ajor 
challenge because it essentially requires a global climatology of cirrus clouds 
that includes microphysical properties. 

There is no measurement system to date that can provide both global 
statistics and microphysical information. Satellite-borne sounders provide 
global coverage, but active spaceborne remote sensing techniques are limited 
in the number of cloud quantities directly measurable, and passive 
instruments have difficulties in retrieving basic parameters such as cloud 
vertical structure, and in the detection of thin cirrus. On the other hand, 
ground-based remote systems such as light radar (lidar) can provide high 
resolution and simultaneous measurements of multiple cloud and 
atmospheric parameters, but with limited coverage. Thus, ground-based 
lidar and satellite measurements provide complementary information, and 
the question we explore here is whether, and how, the knowledge gained 
from lidar measurements can be used to enhance the retrieval of satellite 
data to provide a global, statistical estimate of cirrus cloud radiative forcing. 
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Abstract 

Cirrus measurements obtained with a ground-based polarization Raman 
lidar at 67.9°N in arctic winter reveal a strong correlation between the 
particle optical properties, specifically depolarization ratio and 
extinction-to-backscatter ratio, for ambient cloud temperatures above 
~ — 45°C, and an anti-correlation for colder temperatures. Similar 
correlations are evident in a 2-year midlatitude (53.4°N) cirrus data set. 
Scattering calculations show that the observed dependences can be 
interpreted in terms of the shapes and sizes of the cirrus ice particles. These 
findings suggest a retrieval method for determining cirrus extinction profiles 
from spaceborne lidar polarization data. 
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Three independent cloud optical parameters can be measured with lidar: the 
extinction coefficient, the backscatter coefficient, and the depolarization 
ratio 5 (ratio of elastic backscatter signals with polarization planes 
orthogonal and parallel, respectively, to the linearly polarized lidar radiation 
source). The extinction-to-backscatter (lidar) ratio S and the depolarization 
ratio do not depend on the particle number density, and therefore contain 
range-resolved information about the scattering properties of the cloud 
particles. Simultaneous measurements of lidar and depolarization ratios have 
been obtained with high spectral resolution lidars (4) and polarization 
Raman lidars (5), but correlations between these optical properties and their 
relation to particle microphysical properties have not yet been investigated. 

Figure 1 provides an overview of three cirrus cloud systems measured with 
the polarization Raman lidar of the GKSS Research Center Geesthacht (5) 
above the Swedish research facility Esrange (67.9°N) in January 1997. These 
cases have been selected from cirrus observations on 14 of the 26 days the 
lidar was operated for outstanding data quality and duration of the 
sightings. The evolution of the geometrical structure and S vary strongly 
from case to case. To study the scattering behavior of the cirrus particles 
(<?), height profiles of S and S have been determined for three independent 
observation periods during each cirrus event {7-9). Figure 2 is a plot of 5 
versus 5 for all cirrus altitudes of these 9 measurement segments (A1-A3, 
B1-B3, and C1-C3). The data show two distinct linear correlation patterns. 
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5 and S are negatively correlated when the depolarization ratio is high, and 
they are positively correlated when the depolarization ratio is low. The 
groups merge at S ~ 22 sr with 5 values around 40% (10). The data points of 
individual measurement segments, e. g., B2, are not scattered over both 
subsets of data, but are confined to a relatively small section of one of the 
correlation lines. The correlations strongly suggest that the cirrus 
observations can be given a microphysical interpretation, because S and 5 
are determined by the scattering phase matrix of the cirrus particle ensemble 
which in turn depends on the size and shape of the crystals. We investigate 
this hypothesis below by examining the dependence of the particle properties 
on ambient cloud temperature Tq and also by comparison with calculations. 

Temperature is the principal factor that governs the crystal shape (11), and 
so the observed systematic variation of lidar and depolarization ratios should 
depend on Tq. For this purpose Tq values that are associated with the 
particle optical properties of each data subset have been calculated from 
radiosonde measurements. The results of the analysis are presented in Fig. 3. 
Tq is indeed tightly connected to the optical properties, with Tq being 
coldest for cirrus clouds that are characterized by low S and high 5, whereas 
Tq is warmest for low S and low 5. Above S ~ 22 sr cloud temperatures of 
high-5 and low-5 data groups are not distinguishable. 

Secondly, the observed S and 5 values have been compared with theoretical 
particle optical properties. Because T c < -40°C for the overwhelming 
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majority of the data, the existence of water droplets can be excluded (12), 
and it is assumed that the cirrus clouds consist solely of solid ice particles. 
Figure 2 shows the results of calculations for crystals that have a basically 
column-like or plate-like hexagonal shape (13). In general, theoretical and 
observed optical properties agree well. Column and plate 5 values match the 
values of the high-5 (cold T c ) and low-5 (warm T c ) data subsets, 
particularly for small S. For larger S, the computed optical properties mark 
the boundaries of the measurements. This, and the fact that the absolute 
values of the slopes of the two regression lines are almost identical, imply 
that the observed cirrus clouds consist of mixtures of column- and plate-like 
particles. In natural ice clouds length and aspect ratio (length divided by 
width) of the ice crystals are correlated (14), small crystals tend to have 
aspect ratios close to one, whereas large columns and plates have aspect 
ratios » 1 and < 1, respectively. Applying this correlation to the data in 
Fig. 2, we find that the sizes of the scattering cirrus crystals increase with 
increasing S. This, however, does not necessarily mean that the ice particles 
(6) are larger for large S than for small S : A morphologically complex or 
composite particle, such as an aggregate of crystals of simple spatial 
symmetry, has been shown to scatter similarly to the ensemble of basic 
crystals that formed it (15). So a small ice particle consisting of one or more 
relatively large ice crystals may still have a higher S value than a large ice 
particle formed from relatively small crystals. 


7 


In summary, we propose the following microphysical interpretation of the 
observations. In the coldest regions of the clouds, ice particles with aspect 
ratios ~ 1 form. For increasing T c , the ice particles get larger and consist of 
a mixture of column-like and plate-like crystals, but remain predominantly 
column-like. Around Tc = -45°C, the dominant fundamental shape of the 
crystals gradually shifts from columns to plates. At even warmer 
temperatures in the lower regions of the cirrus clouds, particle size and 
morphological complexity increase with Tq. The high correlation of the low-<5 
data subset may be due to the fact that hexagonal plates and complex 
crystals such as hollow columns, bullet rosettes, and aggregates have almost 
the same scattering signature, if rough crystal surfaces are assumed (16, 17). 
This microphysical interpretation combines in a consistent manner the 
observed cirrus optical properties and temperature data with theoretical 
studies (15, 16), laboratory experiments of ice particle growth (18, 19), in 
situ particle sampling (19, 20), and remote measurements of cirrus clouds at 
other locations (21, 22). 

The previously discussed correlations between the particle optical properties 
of Arctic cirrus clouds, specifically 8 and 5, are an important finding on 
their own, as they have proven to be useful for the interpretation of the 
cirrus data measured over Esrange and for defining constraints to ray-tracing 
calculations of particle scattering for remote sensing applications (13). 
However, the question arises as to whether these same correlations are 
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characteristic of cirrus clouds in general (23). A shortcoming of the Esrange 
measurements is the lack of information about supersaturation, after 7c the 
second principal atmospheric parameter controlling particle growth (11). 
Because the tropospheric absolute humidity was low as confirmed by lidar 
water-vapor measurements (1 7), the supersaturation and particle growth 
rates were probably low also, and the results may not apply to ice particles 
that formed rapidly such as those found in glaciated outflows of convective 
clouds. 

As a first check for the generality of the findings, measurements taken with 
the GKSS Raman lidar (5) at its home base at Geesthacht, Germany 
(53.4°N), between 1994 and 1996 have been analyzed for cirrus cloud 
properties. As compared to the cirrus clouds measured over Esrange, the 
observed cirrus clouds over Geesthacht were optically thinner and appeared 
at higher altitudes. Furthermore, atmospheric conditions at Geesthacht were 
often not favorable for lidar measurements. Therefore only cloud-mean 
values of S and 6 can be determined for each cirrus event. The results are 
shown in Fig. 4. It is evident that 6 and S of the northern midlatitude cirrus 
clouds are correlated as well, and that the functional dependence of S on S is 
the same as for the low-5 data group of the Arctic cirrus clouds. The almost 
complete lack of midlatitude cirrus data with 5 > 40%, even for very cold 
ambient temperatures, suggests differences between the growth conditions of 


the ice particles at the two sites. 
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These promising results should encourage further efforts to even better 
understand the relation between optical and microphysical cirrus properties 
and cirrus-particle formation, and to define the atmospheric conditions 
under which the correlations between depolarization and lidar ratio hold. 
This knowledge would have important implications for generation of an 
advanced global cirrus climatology which would then include information 
about the particle extinction coefficient, one of the most important cloud 
parameters necessary for radiative-transfer calculations in climate studies. 
For experimental reasons this quantity cannot be measured directly with 
spaceborne lidars, but the work presented here suggests an indirect method. 
If an elastic depolarization lidar was deployed on a satellite and if the 5 
versus S correlations were applicable, the depolarization measurement could 
be used to derive an estimate of the height-dependent cirrus lidar ratio 
which in turn would allow' to solve the elastic-backscatter lidar equation for 
the particle extinction coefficient ( 24 )- 
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Figure Captions 

Fig. 1 . Height versus time display of depolarization ratio generated from 
consecutive lidar profiles of cirrus clouds measured over northern Sweden. 
From each cirrus observation three measurement periods (arrows) have been 
chosen for further data analysis. The measurement wavelength is 355 nm. 

Fig. 2. Depolarization ratio versus lidar ratio scatter plot of the cirrus 
measurements Al-3, Bl-3, Cl-3 shown in Fig. 1 (•). The distribution of 
the data points is bifurcated, 5 and S of the two obvious data subsets with 
high 5 (O) and low <5 (O) are strongly correlated (regression lines and 
coefficients in corresponding colors). Calculated 5-versus-S values of 
hexagonal columns (▼) and hexagonal plates (a) suggest that, for small 5, 
cirrus scattering is predominantly by column-like or plate-like particles, 
whereas for larger S scattering contributions of both particle types seem to 
mix. Aspect ratios of columns (1.4, 1.3, 1.5, 2.5, 3.8, 5.9, from left to right) 
and plates (0.3, 0.2, 0.13, 0.08, 0.05, 0.04, from left to right) indicate an 
increase in crystal size w’ith increasing lidar ratio (6). 

Fig. 3. Cloud temperatures of the data subsets with low (red) and high 
depolarization ratios (blue, cf. Fig. 2) averaged over 2.5-sr S bins. Error bars 
indicate statistical uncertainties of the temperature means. Temperature 
profiles are taken from quasi-simultaneous radiosonde ascents at the lidar 
site or at the near-by radiosonde station Lulea. 
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Fig. 4. Depolarization ratio versus lidar ratio scatter plot of cirrus clouds 
measured over Geesthacht, northern Germany (53.4°N). Data represent 
cloud-layer mean values of cirrus events on different days. Regression lines of 
Fig. 2 are shown for comparison. 
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Cirrus measurements obtained with a ground-based polarization Raman lidar at 67.9°N in arctic 
winter reveal a strong correlation between the particle optical properties, specifically 
depolarization ratio and extinction-to-backscatter ratio, for ambient cloud temperatures above 
-45°C, and an anti-correlation for colder temperatures. Similar correlations are evident in a 2- 
year midlatitude (53.4°N) cirrus data set. 

Scattering calculations show that the observed dependences can be related to the microphysical 
properties of the cirrus clouds, specifically to the shapes and sizes of the ice particles. This 
microphysical interpretation is consistent with current theoretical studies and laboratory 
experiments of ice particle growth, with in situ particle sampling, and with remote measurements 
of cirrus clouds at other locations. 

This is the first time correlations between the optical parameters and a relation between the 
optical and microphysical properties of cirrus clouds have been reported. The implications of 
these results are important for atmospheric research: First, the correlations between the particle 
optical properties are useful for defining constraints to ray-tracing calculations of particle 
scattering for remote sensing applications. Second, the observations reveal a new method for 
determining the microphysical properties of cirrus clouds, which can then be used to infer the 
ratio of light transmitted and scattered back to space. This is not only important for similar 
instruments, such as spaceborne lidars which will fly over the next few years, but also for the 
understanding of the general cirrus cloud radiative transfer problem. 


